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Abstract: 

The selfconsistent approach to the two-dimensional Ising model with quenched random 
bonds is extended to the full lattice theory described by four real fermions. A calculation of 
the averaged spin-spin correlation function for large separation of the spins in the disorder 
dominated phase indicates an exponential decay of this quantity and therefore a vanishing 
spontaneous magnetization. The corresponding correlation length is proportional to l/ry 2 , 
where 77 denotes the order parameter of the new phase introduced by Ziegler. 
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0. Introduction 



The 2D-Ising model [1], defined on a square lattice with quenched disorder in the ferro- 
magnetic bonds is defined via the Hamiltonian 

3 <i,j> 

The bond strengths are ferromagnetic random variables < J(i,j) and J(i,j) = Jo- 
< i,j > is a pair of next neighbor sites. 

Whereas it is well confirmed that this system undergoes a phase transition to a ferromag- 
netic ordered state at a lowered critical temperature compared with the pure system, the 
nature of the transition is not completely understood because of the non-applicability of 
the Harris criterion and several conflicting analytical results [2,3,4]. The only exactly solv- 
able model for a disordered ferromagnet is the McCoy- Wu model [5] , where the disorder is 
essentially one-dimensional. For the model with isotropic disorder the famous analysis by 
V.Dotsenko and Vl.S.Dotsenko [4] mapped the system onto the N = O-Gross-Neveu model 
to get the averaged thermodynamic quantities. This is possible because of the free fermion 
representation of the pure model. The authors used the replica trick and a continuum 
field theory, which was analyzed with the momentum space renormalization group (RG-) 
technique. This model is asymptotically free in the infrared, the coupling g (~ disorder 
strength) is marginal irrelevant, the phase diagram is not changed and the specific heat 
diverges at the critical point but more slowly than in the pure model. (~ loglog(|T — T c \~ l ) 
versus log(|T — T c | -1 )). The averaged spin correlation function < SqS u > shows a drastic 
change, namely slow decay: < S S n > ~ exp — |(log logn) 2 . The disorder seems to en- 
hance the correlations instead of destroying them, the critical exponent goes from 1/4 to 
zero. 

An alternative approach uses the bosonized version of the pure model, which makes the spin 
operators local [6,7,8]. The disorder is treated with the same method as in [4]. The results 
confirm the behavior of the specific heat but predict only logarithmic corrections to the 
correlation exponent. An equivalent way is to take advantage of the conformal invariance 
of the pure model [9]. All these methods bear two difficulties: first, the replica trick is 
used and second, they start from the continuum limit associated with the supposed critical 
point. Now, there exist an alternative to the replica trick through a supersymmetric version 
of the effective model [3,10]. This is important because the replica trick is questionable in 
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our case [11,12]. [10] performs a RG-analysis of this model and [3] calculates the saddle 
point structure of the theory after a transformation to composite operators (Q-matrices). 
Both [3] and [10] show finite specific heat in the critical region. Whereas the RG-treatment 
in [10] lacks the correct incorporation of the additive renormalization crucial to extensive 
quantities [13], [3] shows that in the critical region a new saddle point becomes stable 
and governs the thermodynamic behavior of the effective model. This saddle point is 
accompagnied by spontaneous symmetry breaking and a new phase between the ferro- and 
paramagnetic one. The corresponding order parameter stems from a regularization term 
needed for the boson integration in the supersymmetric theory [3,10]. The non- vanishing 
of the order parameter can be proven rigorously [14,15], but there exist two regularizations 
corresponding to different physical situations, one to the Random bond Ising model, the 
other to a system of polymer chains [12]. 

All investigations mentioned so far use a large scale approximation where two of four lattice- 
fermion degrees of freedom are ignored. The RG-approaches even need the continuum limit. 
The difficulty of the continuum limit is connected with the fact that the renormalization 
procedure and the bosonization are not interchangeable. This is the "technical" reason 
for the discrepancy of the findings in [4], respective [5,6,7] regarding the spin correlation 
function [16]. 

Therefore it seems worthwile to extract as much information on the magnetic correlations 
as possible without large scale approximation. This is the aim of the present article. 
Moreover, the technique used avoids both the replica trick and the supersymmetric theory, 
staying close to the original model by direct average of the Greensfunction in the framework 
of a 1/N expansion [2]. The parameter N serves as a bookkeeping device to derive a 
selfconsistent theory which contains the tadpole structure completely (section 1). The 
diagrammatic expansion builds a bridge between the Q-matrix approach and the RG- 
calculations by comparison of the class of Feynman diagrams which are accounted for. 
The tadpoles are usually omitted in field theoretical investigations because they can be 
treated by normal ordering. But this may be dangerous if one does not know the vacuum 
structure, i.e. the phase diagram. With the full 4-fermion lattice theory one can identify 
the diagrams which yield a contribution second order in the coupling to an "exceptional" 
mass term (see section 1), which turns out to be equivalent to the regularization term 
in [3]. Eventually, this term (the order parameter of the new phase) cannot be viewed 
as an artefact of the supersymmetric theory. But the main reason to consider the full 
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theory is the possibility to calculate magnetic correlations which are the subject of the 
many computer simulations performed on the model [17,18,19]. These show very good 
agreement with the thermodynamic predictions of [3] (see [2]; an experimental result for 
the specific heat is also available [20]) 

The disorder average of the square of the spontaneous magnetization in the new phase is 
now calculated in section 2. This is done by determination of the averaged spin correlation 
function for large separation of the spins (section 2). This quantity is related to the 
averaged square of the extensive spontaneous magnetization, i.e. Ai 2 = (^2i(Si)) 2 , which 
is translationally invariant before the disorder average, by means of a Griffith inequality 
[21] and the cluster property [21,22] together with the hypothesis of self-averaging of A4. 
(The spin correlation function itself, whose first moment is calculated in this article needs 
not to be self-averaging [9]. But exponential decay of the magnetic correlations (SiSj) 
is sufficient for the conclusion, that M. vanishes in a fixed sample in the thermodynamic 
limit.) Even if the cluster property is questionable [22], our result (M = 0) is not affected 
by this. However, M 2 should not be confused with the Edwards-Anderson order parameter 
for spin glasses [22]: {Si) 2 , a local quantity, which becomes translationally invariant after 
the average over disorder. Therefore, the present methods do not allow for an estimate 
of this important quantity in the new phase. We use Gaussian disorder to simplify the 
calculations although a bounded distribution of the bond strenght is necessary to keep all 
the bonds ferromagnetic. Nevertheless, the only relevant cummulant in the N = oo limit 
is the second one. In [2] it was shown that the higher cummulants appear only in higher 
order terms in the 1/iV-expansion. 



4 



1. The selfconsistency equation. 

The N-fold replicated partition function for a specific configuration of disorder is 

Z DN = J £>£ expire 
with the euclidean action (Hamiltonian) 1,2 

N 

Hd = £ £ff (r)ff(r) + ft(r)ft(r) + £?(r)tf(r) +ff (r)£(r) + £ 3 a (r)£ 2 Q (r) +&(t)£2(t) 

a, 13 r 

+t ^(r)^(r + e y ) + toft (r)tf (r + e x ) (1) 
(r)ff (r)£?(r + e B ) + <5#V)tf(r)tf(r + e x ). 

The £f are real Grassmann fields defined at the points r of a two dimensional square lattice 
A = Z x Z, e XtV are the unit vectors in the two directions of the lattice, to = tanh Jo/kT 
measures the average bond strength. The disorder variables 5t%P are statistically indepen- 
dent and Gaussian distributed with mean zero and variance = x, y) 



8t^{v)5tf{T') = jjS^SijiSafas + SagSfr]. (2) 

For a given configuration D of disorder one expands the Green's function for a fixed index 



a: 



G'5 D (r,r') = (ef(r)^(r'))A I5 (3) 



With 



one gets 



GS°(r,r') = (4 a (r)^(r'))^=o (4) 



G?f{T,T>)=G?j > + 

oo N m 

E E U.^v k ' Pk M^ k (r k )^(r k + e y ) (5) 

m=l ct k ,/3 k r 1 ...r m 

+St^ft k (r k )^ k (r k + e x )}Z?(r')) 5t=0 . 
Applying the fermionic Wick theorem and averaging over the disorder fields to get the 



averaged Green's function G'(r, r') = G a (r, r') leads to a closed loop of G°-propagators 
with selfcontacts. As was shown previously 2 the only terms contributing to order N° have 
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the tadpole structure. There are essentially four of them (i.e. eight for both axes), depicted 
in fig.l and fig. 2. The left graph of fig.l, for example, corresponds to 



9_ 

N 



N 



E E G ti^ nK'iVi + e y , n)G^(n + e y , r') 

/3 ri 



(6) 



Therefore one deduces a (nonlinear) selfconsistency equation including all tadpoles. Drop- 
ping the replica indices it reads 

G" M (r,r') = G° M (r,rO+ 
9 E E E (- 1 ) 5hjG li^ 'OGjjCri + e(i), n + e^/n + e(j), r') (7) 



ri i • 



with 



and 



1 = 2,2 = 1,3 = 4,4 = 3 

e(l) = -e(2) = -e y 
e(3) = -e(4) = -e x . 

The sign factor in (7) comes from an exchange of Grassmann variables in case of diagonal 
propagators in the loop (see fig2.). Equation (7) determines the "selfenergy" of the effective 
Hamiltonian. In momentum space one gets with H' = G" -1 and H = (G ) -1 



H (p)=H'(p)+gC(p). 



(8) 



G is a 4 x 4 matrix 



/ c 0y 



c 



e lPV Cn, 



with 



CQy 











CQx 
- P -iPx r * 



o \ 





e Wx c x 

CQx ) 



(9) 



-F 

J —7 



dp x dPy 



2tt 



V p "G ,21 (p) 



(10) 



J —IT 



dp x dpy iPx , 
2ix 43(Pj 



Now for x/y-symmetric disorder, cq x = co y and the kernel of the Hamiltonian of the pure 
system reads 1 

/ a 1 1\ 
-a* -1 1 
-110 6 
V -1 -1 -b* 0/ 



#o(*o,p) = ^ 



(11) 



with a = 1 — t e iPy and 6=1 — t e %Px . t = tanh(l//c B T). Therefore, the most general 
ansatz compatible with the selfconsistency equation is 

H'(t,p) = H (t,p) + ^l (12) 

where 1 is the identity matrix. A nonvanishing rj G R prevents H' from getting eigenvalues 
= for any value of p and the divergency of the specific heat in the pure model vanishes. 
In a large scale approximation one recovers just the "externally regularized" model of [3], 
where the possibility of the ?y-term was assumed for different reasons. In fact, a saddle 
point calculation led to a r\ ^ in a narrow region around the (shifted) critical temperature 
which corresponds to t = t c = \[2 — 1. 

In the present consideration of the full lattice theory, the ?y-term is connected with the 
diagonal propagators of fig. 2. These vanish if rj = because the pure Greensfunction 
G%(r,r) = 0. This is clear as we started from a Majorana-field theory with real fermions. 
The ?y-term is forbidden within such a field theory. However, the average over disorder 
gives nonzero contributions to the diagonal propagators second order in the coupling g. 
But the selfconsistency condition allows for rj ^ only in a neighborhood of the critical 
temperature as in [3]. Eventually the model must correspond to an effectice Dime-Held 
theory in this region which allows for the diagonal entries in H'. The breaking of the 
discrete symmetry in [3] can be understood as "spontaneous charge generation" in going 
from real to complex fermions. 

In our approach the condition for a nonvanishing rj deduced from (8) is 
with 

det(r,, t) = det(2#'(r7, t)) = + Ar, 2 + r] 2 (\a\ 2 + \b\ 2 ) - 4Re(a)Re(b) + \a\ 2 \b\ 2 + 4 (14) 
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This expression is minimal for t = t c , therefore it exists a region around the point t = t c 
determined by 

f , 2 rf + bb'+2 

i=a r p te^t) (15) 

where rj ^ 0. This equation, together with 

, , , n f * e^[a*(\b\ 2 + V 2 )-2Re(b)] 

t0 = t + 9 J dp d^J) (16) 

allows to calculate rj(t ,g) and t(t ,g) in this region. The phase boundaries to the "outer 
phase" with rj = are given by 

1= ,/>p ■ W+ \ - , (17) 
i &f(l7 = 0,t(f± g)) 

Equation (16) leads for r\ = to 
to = t + g 



t 3t 2 -l i,l 1 , 3t 2 - 1 

F(-,-;l;/c 2 )-s— — — (18) 



+ 



1 + t 2 4t(l-t 2 )l v 2'2' ' ' *4t(l-t 2 ) 

with F(/c 2 ) the hypergeometric function and k = 4t(l — t 2 )/(l + t 2 ) 2 . Fig. 3 shows the 
effective t as function of to an d the two values £_ and t + . In the next chapter it is shown 
that the spontaneous magnetization vanishes in the new phase characterized by rj ^ 0. 
Therefore, t+ has to be viewed as the transition point from the ferromagnetic to the 
disordered phase. This leads to a lowered critical temperature, as expected. 
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2. The spontaneous magnetization in the disordered phase. 

To get the spontaneous magnetization, we calculate the averaged spin-spin correlation 



function along the x-axis S(n) = (a(O)cx(n)) for large values of the spin distance n measured 
in units of the lattice spacing. This quantity is connected to the magnetization through 
the relation 1 



M 2 = lim S(n) (19) 

n^oo 

The following discussion has two limitations: 

1. We use the effective Hamiltonian derived in section 1., the inverse of the disorder 
averaged Greenfunction for the fermion operators. This leads to the omission of certain 
contributions to the averaged correlation function of the spin operators which are nonlocal 
in the fermions 1 . Nevertheless, this is correct within the N = oo approximation because 
the omitted terms are of higher orders in 1/N. For a specific configuration D of disorder, 
Soin) can be written as ratio of two partition functions 1 

n-l „, 

S D (n)= Ht x (i,0)^ (20) 



i=0 



D 



where in Z' D the bonds along the x-axis between the sites (0, 0) and (n, 0) are modified 
from t x (i,0) to t~ 1 (z,0). This amounts to replace ^H D ^ in the functional integral by 
£H D Z + \m and 

Z l x {l,U) i=o 

Now (exp \£,Q£)z D is expanded 4 in a series of closed loops of products of fermionic prop- 
agators atached to the line between (0,0) and (n, 0) (see Fig. 4). Averaging now these 
quantities leads to contributions of the forms depicted in Figs. 5, 6 and 7. If we replace 
in formulas (20) and (21) t x (i,0) by t and Hp by H' as derived in sect.l, the graphs of 
Fig. 5 and 6. are counted but not those of Fig. 7. However, only the first two are of order 
iV° in the 1/N expansion. This corresponds to the fact that the level of a renormalization 
group treatment is not reached in this consideration on "mean field level", which is suit- 
able to handle the nontrivial phase diagram of the problem. Nevertheless, the graphs of 
fig. 7 can be factorized from the contributions which are accounted for in the expansion of 
the averaged correlation function. Therefore they cannot affect the result of the following 
calculation, i.e. they cannot lead to a nonvanishing spontaneous magnetization in the new 
phase. 
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2. The results are only valid in a small region around t + , i.e. for small rj and welldefined 
Go(t(to: g))- In addition we use Szego's lemma 1 to evaluate the spin correlation function 
of the associated "pure" system, i.e. with the modified coupling t instead of to- This 
gives results only in the limit n — > oo. Nevertheless our result is applicable for large but 
finite n, because in the vicinity of t+ the strong n-dependence of S(n) in the r] ^ region 
is multiplied by an almost constant factor coming from the nonvanishing spontaneous 
magnetization for t < t c . Moreover, the width of the region of validity depends only on 
not on n. The approximations made in the appendix affect the analytical form of the 
dependence of the correlation lenght on rj only (see equation (37)). They have no influence 
on the qualitative behavior of the averaged correlation function. Close to t+ our result is 
correct at least in the sense of the "scaling limit" results obtained in the pure model [1]. 
Following [1] we write for the square of the disorder averaged spin correlation function 



Q' projects on the (3,4)-(i.e. horizontal) sector of the four fermions and is off-diagonal in 
the position space index. 



x(x) means the characteristic function on the intervall [0, n]. Eventually, Q' projects G' on 
a 2(n + l) x 2(n + l) dimensional subspace. Therefore, we have to evaluate the determinant 
of the 2(n + 1) x 2(n + 1) matrix 




(22) 



Qi,j(xi,yi;x 2 ,y2) = Si^8j A S X2 , Xl +ix(xi)S yij0 8 y2fi - 5i A S jy 35 XljX2+1 x(x2)S yi ,o5y 2j o (23) 



tl + (1 - t)G'Q'. 



(24) 








t 



1 



(25) 



1 



t 








Now G' = 



(H (t) +r)l) 1 is expanded to second order in 77, (G = H (t) 1 ), 



G'(t, v )^Go(t)(l- V Go(t)+ V 2 Gl(t)). 



(26) 
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Therefore, the determinant of (24) becomes 



det(C + (1 - t)G (-vG + rfGDQ') 



(27) 



with 

C = tl + (1 - t)G Q' 
Eventually, the determinant factorizes 



S(n) = det Cdet(l + (1 - t)C- l G {-r)G + v G 2 )Q') 



(28) 



(29) 



The matrix C gives the averaged correlation function if rj would be absent. Because we 
are in the region where the effective coupling t is greater than t c , det C reaches a constant 
nonzero value in the limit n — ► oo: 



lim det C = [1 

n^oo 



(sinh2/3)- 4 ] 1 / 2 



(30). 



(3 = atanh(t). Equation (30) gives the averaged square of the magnetization in the "outer" 
phase, i.e. for t < t- or t > t+. It leads to a nonvanishing magnetization for t > t + 
as in the pure system, albeit the functional dependence of the effective t on to has to 
be taken into account. However, in the framework of the iV — > oo limit the averaged 
magnetization exhibits a finite jump at the point t + . That means, the phase transition from 
the magnetically ordered to the disordered state occurs at a temperature corresponding to 
t+. The spontaneous magnetization vanishes for t < t+, whereas it is different from zero 
for t > t + . To see this, we calculate the second determinant in (29) for < rj <^ 1. 
C has the following block structure 1 



/ 






1 



\ 







V 



i 






J n-1 



At 



(31) 



A is a n x n Toeplitz matrix with the elements 

" n dp 



A 



2n 



f(p)e 



ip(j-k) 



(32) 
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where we choose the indices j, k to run from —(n — l)/2 to (n — l)/2, and 

bj = |* ^f( P )e^ n -^ (33) 

for j G [0,n — l]. /(p) is an unimodular function, parametrized by t and its dual t, 
t=(l-t)/(l + t): 

/(p) = ( (te^-t)(tte^-l) J • (34) 

The second factor in equation (28) can now be written in the form 

exp(Tr log(l + (1 - ^C^GoC-t/Go + rj 2 G 2 )Q')) (35) 

The logarithm can be expanded into powers of rj, because we are not in the critical regime 
and therefore Go(t) is a bounded operator in its domain of definition and has finite operator 
norm. By an appropriate choice of rj, (close enough to zero), the expansion is welldefined. 
In first approximation higher terms then quadratic are neglected. In the appendix it is 
shown that the term proportional to rj 

Tril-^C^GlQ' (36) 

vanishes, and the term proportional to rj 2 is negative definite, yielding finally for large n 



S(n) ~ const, x exp[— (rj 7(£))n] (37) 

where ^(t) is a positive constant, depending on to and the disorder strenght g. One sees 
that the disorder averaged spin correlation function decays exponentially for an arbitrary 
nonzero value of rj. Eventually, the averaged spontaneous magnetization vanishes in the 
phase with rj ^ 0. 
Conclusions 

In this article, the phase transition from the ferromagnetic ordered phase to the disor- 
der dominated phase in the two dimensional Random Bond Ising Model is investigated. 
Whereas in previous approaches 2 ' 3 ' 4 a large scale approximation in the vicinity of an as- 
sumed critical region allowed to reduce the number of independent fermion degrees of 
freedom from four to two, we treat the full lattice theory in a selfconsistent way in the 
framework of a 1/N expansion. This expansion is not supposed to correspond to some spe- 
cial physical properties of the system, but serves as a bookkeeping device, which allows to 
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extract the most important "mean field" contributions usually omitted in renormalization 
group approaches (i.e. tadpoles). The most general selfconsistent ansatz for the effective 
Hamiltonian has the same structure as the saddle point solution of the Q- matrix theory 3 , 
although the reason for introducing a possibly nonvanishing r\ is different: In the Q-matrix 
theory one uses supersymmetric fields and the ?y-term regularizes the integration over the 
bosons. After averaging over the disorder, r\ becomes nonzero in a certain region due to 
spontaneous symmetry breaking. The order parameter in the symmetry broken phase can 
therefore be identified with rj, leading to a continous phase transition from the ferromag- 
netic to the disorder dominated phase as well as from this phase to the paramagnetic one. 
The question arises, whether the spontanous magnetization in the new phase vanishes or 
not. To this end, we calculated the asymptotic, (n — > oo), behavior of the disorder av- 
eraged spin correlation function in a region around the transition line from the magnetic 
ordered phase to the new phase. Whereas the magnetization on the ferromagnetic side is 
finite corresponding to T < T c , the correlation function decays exponentially 

S(n) ~ exp — n/£ 

with 

4 7(*)<7 2 

That means, rf plays the role of an inverse correlation length. However, it is not clear 
what kind of operators becomes massless as £ goes to infinity. The Q-matrix analysis 
shows a massless mode of the local composite Q-operators at t + . Eventually, a rigorous 
RG-treatment of the transition point has to deal with these objects. The temperature 
dependence of the magnetization in the ferromagnetic phase can be calculated by means 
of equations (18) and (30). The appearance of ?y 2 and not rj in the expression for the 
averaged correlation function is due to the fact, that the selfconsistent equation allows two 
real solutions with rj positive or negative, whereas the physics must not depend on the sign 
of r\. This corresponds to the freedom to choose the sign of the regularization term in the 
supersymmetric theory. If one relates the averaged asymptotic correlation function to the 
averaged square of the spontaneous magnetization in the usual way, it follows, that the 
new phase shows no long range magnetic order. Therefore, the only possibility to discern 
it from the paramagnetic region consists in the investigation of the relaxational dynamics 
for the spins, which will be discussed elsewhere. However, the present article says nothing 
about the transition from the new phase to the paramagnetic one. It is even possible, 
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that something happens in the new phase itself, when the effective t reaches t c . In the 
framework of the N — > oo, resp. saddlepoint calculation, the spontaneous magnetization 
exhibits a finite jump at the transition from the ordered to the new phase. Whether 
this remains correct if one takes into account fluctuations can perhaps be answered via a 
renormalization group treatment of the theory on the ferromagnetic side. 
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Appendix 

To calculate the term proportional to r\ in equation (36), we write Go(p) as 



G = 



Oil 


Pi 


e 


5 \ 


-PI 


a\ 


-5* 


e* 


-e* 


5 




P2 


-5* 


— e 


-PI 


a* 2 J 



with 



detH 



«i = b — b* 
ct 2 = a* — a 
f3 x = (b + b*)-abb* 
02 = (a + a*) — baa* 
e = ab* - 2 
S = ab - 2 



(Al) 



(A2) 



with the conventions of equation (11). The (3,4)- (i.e. horizontal) sector of Gq(p) rea ds 
then 

J_^-(|a 2 | 2 + |/3 2 | 2 + |e| 2 + |5| 2 ) 

-(|a 2 | 2 + |/3 2 | 2 + |e| 2 + |5| 2 ) 



detif 2 







(A3) 



Therefore, GqQ' has the form 



* 
* 



(AA) 



where and * are (n + 1) x (n + 1) matrices, multiplying from the left with C _1 yields a 
matrix with vanishing diagonal entries. It follows that the contribution proportional to rj 
vanishes. 

To calculate the quadratic term, we write 



TriC^GlQ') = Tr(G 2 (C - fl)C -1 ) = fr(G 2 ) - f r(G 2 flC'- 1 ) 



(A5) 



and Tr means the trace in the 2(n + l)x2(n + l) dimensional subspace of the horizontal 
sector under consideration. To give an estimate for the matrix C -1 , we use the fact that 
the matrix A in equation (31) is "almost unitary", in the sense that 



AA^ = 1 + -B 



n 



(A6) 
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where B has bounded matrix elements. This is due to the fact that 

(n-l)/2 

E 9(p)e ipj = 9(0) + 0(l/n) 

j = -(n-l)/2 



(A7) 



for a bounded function g(p), and /(p) in equation (32) is of unit modulus. Then it is 
possible to write for C _1 



C 



-l 



At 



d 






dn-l 1 



\ 



1 <-l 





n 



(A8) 



with 



With the definition 



n-l 



d J = ~Yl 6 *4_(n-l)/2,j-(»-l)/2 = °(V«) 



049) 



i=0 



ip(j-k) 



for j, fee [0, n], we get 

Tr^fiC- 1 ) = 2G^ nn + 2ntz jT ^-h(p)Im(f(p)) + 0(1) 



(A10) 



(All) 



The contribution from the integral leads to at most an oscillating factor for the averaged 
spin correlation function, but actually it vanishes due to the symmetry properties of h(p) 
and f(p). It follows that the relevant contribution is the first term in equation (A5), 
namely 

frG 2 = 2(n + l)G ,nn. (A12) 
Finally one gets for the (positive definite) 7(£) of equation (37) 



7(0=2(1-*) 



dp x dp y \a 2 \ 2 + \(3 2 \ 2 + |e| 2 + \5\ 



det H (t) 2 



(A13) 
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Figure Captions: 

Fig. 1: Two off diagonal contributions to the propagator, coming from averaging over 
bond disorder in y-direction. The left graph corresponds to equation (6). a and (3 are 
indices of the N colors, a is fixed whereas (3 runs from 1 to N. 
Fig. 2: Two diagonal contributions. These vanish if r\ = 0. 

Fig. 3: A two-propagator term in the expansion of (exp \£Qti)z D - The horizontal line 
indicates the section of the x-axis between the two spins. The distance between them in 
units of the lattice spacing is n. 

Fig. 4: A term in the average of the expression (exp ^Q0z D , taken into account using 
the effective G' from section 1. a is a fixed replica index and (3 runs from 1 to N . 
Fig. 5: A term taken into account by replacing to by t. (3 and 7 run from 1 to N . 
Fig. 6: A term which is omitted in the N — > 00 approximation. 

Fig. 7: The effective bond strenght t as function of t . The vertical lines denote t_ = 0.18, 
resp. t + = 0.71, the new phase lies between them. The value of g is 0.3. 
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